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Abstract. The accretion process onto a magnetized isolated neutron star, which captures material from the inter- 
stellar medium, is discussed. The evolutionary track of such a star can be presented as a sequence of four states: 
ejector, supersonic propeller, subsonic propeller, and steady accretor. I show that subsonic propeller — > accretor 
transition does not occur as long as the magnetic field of the star is strong enough to control the accretion flow 
in the stellar vicinity. During the subsonic propeller state the accretion rate onto the stellar surface is limited 
to the rate of plasma diffusion into its magnetosphere. The diffusion rate is at least three orders of magnitude 
smaller than the maximum possible mass capture rate by the star. Therefore, the expected accretion luminosity 
of magnetized isolated neutron stars is at least three orders of magnitude smaller than that previously evaluated. 
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1. Introduction 

The observational appearance of isolated neutron stars 
(INSs), which accrete material from the interstellar 
medium, is a subject of intensive theoretical investiga- 
tion (e.g. Shvartsman [Tn7Dl Treves & Colpi HWR Treves 
et al. I2HDDI Toropina et al. I2UUT1 Popov et al. I2000ajl . 
It is presently established that the initially fast rotating 
INSs, under certain conditions, are able to switch their 
state from ejector to propeller within the Hubble time 
t < t-R « 10 10 yr, and even to reach the periods at which 
the corotational radius, 

R cm = (GM ns /cu 2 ) 1/3 , 

exceeds the magnetospheric radius of the star, 

Rm = (M 2 MW2GA/ ns ) 2/7 . (1) 

Here M ns , lu — 2tt/P s and p are the mass, the angular 
velocity, and the dipole magnetic moment of the neutron 
star, respectively. M c is the mass capture rate by the star 
from the interstellar medium, which can be expressed as 

ttRIrVoo. (2) 



Limiting p < 10 gem - ' 3 and Vx, > V s , where V s 
is the thermal velocity in the interstellar material, Popov 
et al. (2000a) have estimated the maximum possible rate 
of mass capture by an isolated neutron star as 



M c <, M n 



10 12 p. 



24 m 2 V 6 gs" 



(4) 



M, 



Voo is the relative velocity of the neutron star to the sur- 
rounding material, p is the density of the interstellar ma- 
terial, and R a is the accretion radius of the star, 



R a = 2GM ns /V c 



(3) 
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Here m is the mass of the star expressed in units of Mq, 
p_24 = /5/10- 24 gcm- 3 and V 6 = cms" 1 . Then, 

assuming that all material captured by the star is accreted 
onto its surface, they have suggested that the accretion 
luminosity of INSs can be as high as M max GM ns / R ns ~ 
10 32 ergs _1 . 

The analysis of the validity of this assumption is the 
subject of the present paper. I show that under the con- 
ditions of interest the mass accretion rate onto the stellar 
surface is limited by the rate of plasma penetration into 
the stellar magnetosphere, which is significantly smaller 
than Af max . Therefore, the maximum possible accretion 
luminosity of a magnetized isolated neutron star proves to 
be limited to a few xl0 28 ergs _1 . The evolutionary tracks 
of isolated neutron stars are briefly reviewed in the next 
section. The rate of plasma penetration into the mag- 
netosphere and, correspondingly, the mass accretion rate 
onto the surface of isolated neutron stars are evaluated in 
Section 3. Here, I also show that the rate of mass accre- 
tion onto the surface of an isolated neutron star remains 
smaller than M c as long as R m 3> R n s- The main conclu- 
sions are summarized in Section 4. 
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2. Evolutionary tracks of INSs 



As was first recognized by Shvartsman 1|1970[1 . the evolu- 
tionary track of a rotating magnetized neutron star can 
be presented in the form of the following sequence of 
its states: ejector — > propeller — > accretor. Within this 
scheme, the rotational rate of a newly born fast rotating 
neutron star decreases, initially by the generation of the 
magneto-dipole waves and ejection of relativistic particles 
(pulsar-like spin-down), and later by means of the interac- 
tion between its magnetosphere and the surrounding ma- 
terial (propeller spin-down) . The first state transition oc- 
curs when the pressure of the material ejected by the star 
can no longer balance the pressure of the surrounding gas, 
and the latter, penetrating into the accretion radius of the 
star, interacts with the stellar magnetosphere. A detailed 
analysis of this state transition with respect to INSs was 
presented by Popov et al. {2000b). 

The spin evolution of a spherically accreting strongly 
magnetized neutron star in the state of propeller has 
been investigated by Davies, Fabian & Pringle (1979) and 
Davies & Pringle ljl981|l . As they shown, two sub-states 
of the propeller state can be distinguished: the supersonic 
and subsonic propeller. In both cases the neutron star is 
spinning down due to the interaction between its magneto- 
sphere and the surrounding gas. As a result of this interac- 
tion, the star's magnetosphere is surrounded by a spherical 
quasi-static atmosphere, in which the plasma temperature 
is of the order of the free-fall temperature 



T ff (r) = (GM ns m p )/(kr) 



(5) 



Here G, m p , and k are the gravitational constant, the 
proton mass, and the Boltzmann constant, respectively. 

The atmosphere is extended from the magnetospheric 
boundary up to the accretion radius of the neutron star. 
The rotational energy loss by the neutron star is convected 
up through the atmosphere by the turbulent motion and 
lost through its outer boundary. 

The formation of the atmosphere in the first approxi- 
mation prevents the surrounding gas from penetrating to 
within the accretion radius of the star. As the neutron star 
moves through the interstellar medium, the interstellar gas 
overflows the outer edge of the atmosphere with the rate 
M c (see Eq.|21), which is traditionally called the strength of 
the stellar wind and denotes the maximum possible mass 
capture rate by the neutron star. 

2.1. Supersonic propeller 

As long as the angular velocity of the neutron star is large 
enough for the corotational radius to be smaller than the 
magnetospheric radius, the star is in the centrifugal inhi- 
bition regime (i.e. the centrifugal acceleration at the mag- 
netospheric boundary, us 2 R m , dominates the gravitational 
acceleration, GM ns / R^). The centrifugal inhibition is not 
effective only within the bases of the corotational cylinder. 
However, the accretion of material onto the stellar surface 
through these regions does occur only if the the angle be- 



tween the magnetic and rotational axes is small enough 
(see Ikhsanov 2001c) and if the magnetic field of the star 
is weak enough for the magnetospheric radius to exceed 
the stellar radius only by a factor of 2-3 (for discussion see 
Toropin et al. 119991 and Romanova et al. l2002|l . Otherwise, 
the accretion power is significantly smaller than the spin- 
down power due to propeller action by the fast rotating 
star. 

Except the bases of the corotational cylinder, the linear 
velocity at the boundary of the magnetosphere, which is 
co-rotating with the star, in this case is larger than the 
sound speed in the atmospheric plasma. That is why this 
state is usually refereed to as a supersonic propeller (see 
also Ikhsanov 2002). 

2.2. Subsonic propeller 

As the star is spinning down, its corotational radius in- 
creases and reaches the magnetospheric radius when P s = 
P c d, where 



P,, 



23 ffl m-^ 



(6) 



Here M15 is the strength of the stellar wind expressed in 
units of 10 15 gs _1 , and is the dipole magnetic moment 
of a neutron star expressed in units of 10 30 Gem 3 . 

Under the condition P s > P c d the centrifugal barrier is 
not effective: the atmospheric plasma, penetrating into the 
magnetic field of the star, is able to flow along the mag- 
netic field lines and to accrete onto the stellar surface. 
However, as shown by Arons & Lea l|197fi[) and Eisner 
& Lamb l)197ri|l . the rate of plasma penetration into the 
magnetosphere of a spherically accreting strongly mag- 
netized neutron star can be as high as M c only if the 
magnetospheric boundary is unstable with respect to in- 
terchange (e.g. Rayleigh- Taylor) instabilities. Otherwise, 
the rate of plasma penetration is limited to the diffusion 
rate, which is a few orders of magnitude smaller than M c 
(see Eq. I10|l . For instability to occur the sign of the ef- 
fective gravitational acceleration at the magnetospheric 
boundary should be positive: 



5cff 



GM ns 
RIM 



cosf 



XRm) 



Rcurv(9) 



> 0. 



(7) 



Here i? C urv is the curvature radius of the field lines, 9 is 
the angle between the radius vector and the outward nor- 
mal to the magnetospheric boundary and Vr^-Rm) is the 
ion thermal velocity of the accreting plasma at the bound- 
ary. For the case of the equilibrium magnetospheric shape 
derived by Arons & Lea l|1976f) . the condition Q can be 
expressed in terms of the plasma temperature at the mag- 
netospheric boundary as 



T <T„ 



0.3 T ff . 



(8) 



This indicates that the condition i? m < R cor is necessary, 
but not sufficient for the effective accretion (with the rate 
of ~ M c ) onto the surface of neutron star to start. In 
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addition, it is required that the cooling of plasma at the 
base of the atmosphere is more effective than the heating. 

As shown by Davies & Pringle 1)1981[) . the cooling of 
the atmospheric plasma is governed by the bremsstrahlung 
radiation and the convective motion. For these processes 
to dominate the energy input into the atmosphere due to 
the propeller action by the star, the spin period of the star 
should be P s <; f^ r , where Pb r is a so-called break period, 
which according to Ikhsanov ( 20 01af) is 

P br ~ 450 M ^ /21 Mf 5 5/7 m- 4 / 21 s. (9) 

Under the conditions of interest, the break period signif- 
icantly exceeds P c d- This means that the state transition 
supersonic propeller — > steady accretor may occur only via 
an additional intermediate state, which is called the sub- 
sonic propeller. This term reflects the fact that the ro- 
tation velocity of the magnetosphere during this stage is 
smaller than the thermal velocity in the surrounding gas. 

If the propeller action were the only source of heating 
of the atmospheric plasma, the magnetospheric boundary 
of the neutron star would be able to switch its state from 
subsonic propeller to accretor as its spin period reaches 
Pbr- However, as shown below, an additional heating of 
the atmospheric plasma occurs due to a radial plasma drift 
through the atmosphere. This additional heating mecha- 
nism turns out to be not effective if a star is situated in a 
relatively strong stellar wind, but in the opposite case it 
must be taken into account. 

3. Diffusion-driven accretion 

Although the interchange instabilities of the magneto- 
spheric boundary during the subsonic propeller state are 
suppressed, the 'magnetic gates' are not closed completely: 
the atmospheric plasma is able to penetrate into the stel- 
lar magnetic field due to diffusion. The diffusion rate is 
limited to (see e.g. Ikhsanov 12001 bf) 

m b <; lo^gs- Co 1 . 7 : 2 M3-o 1/14 - 1/7 (loi^r) 11/14 -( 10 ) 

Here (o.i = C/0-1 is the diffusion efficiency, which is nor- 
malized following the results of experiments on the nuclear 
fusion (e.g. Hamasaki et al. 11974(1 and the measurements 
of the solar wind penetrating the magnetosphere of the 
Earth (Gosling et al. lTMT|> . 

This means that the hot atmosphere surrounding the 
magnetosphere of the star in the subsonic propeller state 
cannot be purely stationary. For the atmosphere to remain 
in an equilibrium state, the amount of material flowing 
from its base into the magnetosphere must be compen- 
sated for by the same amount of material coming into the 
atmosphere through its outer boundary. Thus the radial 
drift of plasma through the atmosphere with the rate Mb 
and the velocity 

V dr = (M B /M C ) V s (11) 
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towards the neutron star is expected. 

The structure of the atmosphere with the radial 
plasma drift can be explain in terms of the quasi- 
stationary model of Davies & Pringle <|1981[1 as long as 
the characteristic time of the accretion process, 

*dr - R/V dr , (12) 

is larger than the characteristic time of turbulent motions 
t t = R/V%. According to Davies & Pringle l|1981(l . the 
velocity of turbulence motions can be approximated as 
Vt ~ (P m /P) 1/6 wP m . Therefore, the condition i t < im- 
proves be satisfied if the spin period of the star is P s <C 
P qs , where 

P q8 , B * 1.7 x 10 5 CoT! 72 MT 5 3/14 m~ 6/r a4o /14 s. (13) 

Thus, as long as the spin period of the star is P c d ^ P s <C 
Pqs.B, the magnetosphere of the star is surrounded by a 
hot (T ~ Tff ) atmosphere in which the plasma pressure is 
p oc P~ 5 / 2 , the sound speed is c s oc P" 1 / 2 and the number 
density is n oc P -3 / 2 (for discussion see Davies & Pringle 

rrwrji . 

The radial drift of plasma through the atmosphere to- 
wards the neutron star leads to the release of the accretion 
(potential) energy, which is mainly spent in heating the at- 
mospheric plasma. The heating rate due to this process is 
L dl (R) ~ M B GM ns / P. This value is small in comparison 
with that of the spin-down power (see Davies & Pringle 

rn%n E q . 3.2.4), 

L ssp = 1.2 x 10 36 M30 m^PT 3 ergs" 1 , (14) 
as long as the spin period of the star is P s 55 P aC c,B, where 

P aC c,B - 450 Co 1 . 7 ! 6 M30 /42 S. (15) 

However, under the condition P s > P acc the heating of the 
atmospheric plasma is governed by the accretion power. 

The characteristic time of the heating due to accretion 
process is tdr- On the other hand, the cooling time of the 
atmospheric plasma due to the bremsstrahlung emission 
is 

t cool « i br ~ 633 s ( J^) (l^— 3 ) 1 • 

where n is the number density of the atmospheric plasma, 
which can be expressed as n(R) = n(P m )(P m /P) 3 / 2 , 
where 

= ^Sk T (R m y (17) 

Hence, for the temperature of the atmospheric plasma to 
be smaller than the critical value, T cr , the following con- 
dition should be satisfied 

icool(Pa) < tdi(Ra). (18) 

I require this condition to be satisfied at the outer bound- 
ary of the atmosphere since t c oc P and t& x oc P 3 / 2 . 
Combining Eqs. {TTJ, (JEJ, (JTHJ), and iTfTjl. one finds that 
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the bremsstrahlung cooling dominates the heating due to 
accretion power only if the strength of the stellar wind is 
M c ^ M , where 



M ~ 3 x 10 



14 



CoT^—Z-^gS-. 



This means that a magnetized isolated neutron star is able 
to switch its state from subsonic propeller to steady accre- 
tor only if the mass capture rate by this star from the 
interstellar medium is M c <; M . Otherwise, the corre- 
sponding state transition does not occur and the star re- 
mains surrounded by the hot atmosphere. In this case the 
mass accretion rate onto the stellar surface is limited to 
M ^ Mb (see Eq. II 01) and, correspondingly, the accretion 
luminosity is L x < £ m ax, where 



R7 L M- 



11/14 



ergs 



(20) 



Here i?g is the radius of the neutron star expressed in units 
of 10 6 cm, and M u = M c /10 14 g s" 1 . 



4. Discussion 

Application of our findings to the case of INSs leads to 
the following conclusions. First, comparing Eqs. Q and 
<|19|) one finds that under the conditions of interest the 
maximum possible strength of the stellar wind of an iso- 
lated neutron star, M max , is smaller (at least by a factor 
of 3) than Mq. Therefore, the interchange instabilities of 
the magnetospheric boundary of these stars are suppressed 
and the plasma penetration from the base of the hot at- 
mosphere into the stellar magnetic field is governed by 
the Bohm diffusion. In this case the accretion luminosity 
of INSs, whose age is smaller than the characteristic time 
of the magnetic field decay, t m fd, is limited to 



L x {t < t mfd ) = Mi 



GM n 

Rns 



3 x 10 28 ergs 



(21) 



>l/2 -1/14 
''O.l ^30 



11/14 



lO^gs" 1 



This radiation is to be emitted from local regions situated 
at the magnetic poles of the star. That is why the star 
could be observed as a low luminous pulsing X-ray source. 
The minimum period of this pulsar is (see Eq. EJ) 



Pn 



450 /LtgQ 7 m" 



5 / 7 Mi 



-3/7 



(22) 



Second, as the age of an INS becomes comparable with 
t m fd, the magnetic field is almost unable to control the ac- 
cretion flow in the stellar vicinity and the direct accretion 
onto the stellar surface occurs. The accretion luminosity 
of the star in this case is L x = M c GM ns / R ns , i.e. by a 
factor of M c /Mb larger than the accretion luminosity of 
a neutron star in the state of the subsonic propeller. This 
radiation is to be emitted from the whole surface of the 
star and, therefore, the pulsations are not expected to be 
observed. 



Finally, for M c < 10 12 gs , the spin-down time scale 
of a neutron star in the subsonic propeller state is (see 
Davies & Pringle ITMn Eq. 3.2.5) 



(19) t a 



> 



1U C, 0A m fi 30 M 12 r e yr, 



(23) 



where Pq = P a cc,B/10 6 s (see Eq. I15|) . This time scale is 
comparable with the characteristic time scale of the mag- 
netic field decay (see e.g. Urpin et al. H99fi|l . Hence, the hot 
atmosphere surrounding the neutron star in the propeller 
state can be treated within the quasi-static approximation 
suggested by Davies & Pringle (1981). 
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